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Abstract 

In a 2+1-dimensional pure LGT at finite temperature the critical coupling for the 
deconfinement transition scales as (5c{nt) = Jc^t + ui, where n* is the number of 
links in the “time-like” direction of the symmetric lattice. We study the effective 
action for the Polyakov loop obtained by neglecting the space-like plaquettes, and 
we are able to compute analytically in this context the coefficient ai for any SU(A) 
gauge group; the value of Jc is instead obtained from the effective action by means of 
(improved) mean field techniques. Both coefficients have already been calculated in 
the large N limit in a previous paper. The results are in very good agreement with 
the existing Monte Carlo simulations. This fact supports the conjecture that, in the 
2+1-dimensional theory, space-like plaquettes have little influence on the dynamics 
of the Polyakov loops in the deconfined phase. 


1 Introduction and setting 

Let us consider a finite temperature lattice gauge theory (LGT) with gauge group 
SU(A), defined on a 2 + 1 dimensional cubic lattice. The “time” dimension is 
compactified with periodic boundary conditions and compactification length T = 
NtOt, where A/ is the number of lattice spacings in the time direction, and at 
the corresponding lattice spacing; T is the temperature. We consider in general 
asymmetric lattices; denoting by the lattice spacing in the space directions, the 
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( 1 ) 


asymmetry parameter p is defined as: 

p=^ = TNtas . 

at 

Moreover we introduce different bare couplings jSt and j3s in the time and space 
directions. The Wilson action is thenQ 

Sw = NY^ [^AReTrGo*(5^,t)+/5.ReTrG'i2(f,t)) , (2) 

X^t \2=1 / 

where Got and Gu are the time-like and space-like plaquette variables: 

Goi{x,t) = V{x,t)Ui{x,t + l)V^x+ i,t)U-{x,t) , 

Gi 2 {x,t) = Ui{x,t)U 2 {x+ i,t)Ul{x+ ],t)Ul{x,t) . (3) 

We have denoted by Ui{x, t) the link variables in the space directions and by V{x, t) 
the link variables in the time direction. 

For a given p the relation between the couplings jSt and Ps can be obtained by 
requiring that the Wilson action reproduces in the continuum limit the usual 
Yang-Mills action. At the classical level, this leads to the relations 

i^ = “*vTA. p = (4) 

which show that in d = 2 the coupling constant g'^ has the dimensions of a mass. 
This is the main difference with the d = 3 case and it is the reason for which the 
d = 2 models are much easier to study. In fact in this case the coupling constant 
itself sets the overall mass scale. Near the continuum limit a physical observable, like 
the critical temperature T^, with the dimensions of a mass can be written, according 
to the renormalization group equations, as a series in powers of g^. Hence in terms 
of the coupling^ j3t we have 

atTc = ^ ^ H-; (5) 

Pt Pt 

with Cl and C 2 suitable constants. The critical temperature is obtained by looking 
at the critical coupling /3c at which the deconhnement transition occurs. Since the 
lattice size in the time direction is W = we can rephrase eq. (||) as a scaling 
law for the critical coupling (3^ as a function of Np. 

_= JciVt + ai... , (6) 

^The notations throughout this paper are the same as the ones of Ref. Q, to which we refer 
for a more detailed exposition of the general setting of the theory. In particular we choose the 
normalization of the couplings (3s and (St so as to have a smooth large N limit, since we will 
compare the results at finite N with the ones obtained in the large N limit. 

^It is not mandatory to measure in units of at and to choose Pt instead of Ps in eq. (H), but 
these choices simplify the comparison with the results of Sec. 
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where we have introduced the two constants Jc = ci and oi = ^ for future conve¬ 
nience. Finding the precise values of these two constants is one of the aims of the 
present paper. 

It is clear from eqs. i) and (^) that equivalent regularizations with different 
values of p require different values of iV*. Hence, to maintain the equivalence, Nt 
must be a function of p according to eq. ©• 

Among all these equivalent regularizations a particular role is played by the 
symmetric one, which is dehned by 




2 

oiV 


(7) 


(from now on we shall distinguish the symmetric regularization from the asymmetric 
ones by eliminating the subscripts t and s m. {3 and a). By comparing eq.s (|^,^,|^) 
we see that all regularizations are equivalent to the symmetric one provided 


p = l>p, = -, 

(8) 

p 


Nt{p) = pnt , 

(9) 


where rit is the number of links in the time direction with a symmetric regularization; 
nt = Nt{p=l). 

Notice however that these equivalence relations have been derived in the naive or 
“classical” continuum limit, and quantum corrections are in general present. They 
lead in the large-p limit to the following expressions: 

(3t = p{f3 + C^r) + Or + • • • ) = — - - H—f + ... . (10) 

p p^ 

The constants encode the quantum corrections to eq. (||). These corrections 
were studied and calculated in the (3-1-1) dimensional case by F. Karsch in |^. 

It is well known that, as a consequence of the periodic boundary conditions in 
the time direction, the Wilson action enjoys a new Zjv global symmetry. Moreover, 
it becomes possible to dehne gauge invariant observables which are topologically 
non-trivial, such as the Polyakov loop| 


Nt 

P{x) = TrP(x) = Tr n V{x, t) . (11) 

t=i 

The Polyakov loop is a natural order parameter for the Zat symmetry. In the 
high temperature, deconfined phase (T > Tc) we have {P{x)) ^ 0, while in the 
low temperature, confining domain (T < Tc), the symmetry is unbroken, namely 

(-p(n> - 0- _ 

^We will refer to the untraced quantity P{x) as the Polyakov line. 
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The way to approach the description of the deconfinement transition is to con¬ 
struct an effective action for the relevant degrees of freedom, namely the Polyakov 
loops. This would require to perform the exact integration over all the gauge de¬ 
grees of freedom on the space-like links, a clearly impossible task. In H, I we 
discussed a perturbative approach, in which the effective action for the Polyakov 
loop is calculated exactly order by order in an expansion in powers of j3s. Except 
for the theory with N = 2 { see P|) only the zeroth order of this expansion, which 
corresponds to neglecting the space-like plaquettes altogether, has been studied and 
used to determine the critical temperature. The point is that the approximation 
of throwing away the space-like plaquettes appears to be a rather crude one if one 
wants to go beyond purely qualitative results. It is true however that far away from 
the continuum limit, that is for very small Nt, the deconfinement transition occurs 
for small values of (3 (we refer here to a symmetric lattice) and that the couplings 
among Polyakov loops induced by the space-like plaquettes, beeing of higher order 
in /3, can be safely neglected. As we proceed towards the continuum limit (large 
(3) there is no obvious reason why the space-like plaquettes should be uninfluential. 
In fact we know that in 3 -I- 1 dimensions the effect of the space-like plaquettes is 
crucial, as the rescaling of the critical coupling with Nt, which is linear in their ab¬ 
sence, becomes logarithmic, according to asymptotic freedom, if they are taken into 
account. The situation is quite different in 2 -|- 1 dimensions, where the rescaling 
of the critical coupling is linear both with and without space-like plaquettes. This 
is clearly related to the fact that the dimensional coupling constant sets the mass 
scale of the theory. So in 2 -|- 1 dimensions we are in a situation where the zeroth 
order approximation is on one hand a good approximation of the full theory for very 
small Nt, and on the other hand it has the same scaling properties of the full theory 
near the continuum limit. It is natural to conjecture at this point that the effective 
action for the Polyakov loop obtained by neglecting the space-like plaquettes is a 
good approximation even in the continuum limit and up to the temperature where 
the deconhnement transition occurs. Checking this conjecture by comparing the 
analytic predictions based on the zeroth order approximation with the Monte Carlo 
simulations based on the full theory is the purpose of this paper. 

This paper is organized as follows. In the next section we shall construct the 
effective action for the Polyakov loop mentioned above. This section contains ma¬ 
terial which has been already presented elsewhere |]l|, and has been inserted so as 
to make this paper self-contained. In Sec. 3 we shall obtain in the framework of 
the /3s = 0 approximation the scaling law (^ and the explicit value of Oi coefficient. 
The exact calculation of ai is the main new analytic result of the paper. In the 
last section we shall estimate by means of suitably improved mean held techniques 
the Jc coefficient and compare our results with those obtained with Monte Carlo 
simulations. The agreement of the analytic predictions with the Monte Carlo sim¬ 
ulations, which seems to conhrm the above conjecture, is the main physical result 
of our analysis. 
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2 The effective action 


In order to derive the effective action for the Polyakov loops at /3s = 0, one starts 
from the time-like part of the Wilson action (j^) and performs the character expan- 
sionQ 

e5w(/3.=0) ^ \^Y.drDrmXr{Goi{x,t))\ . (12) 

x^t^i Ip J 

In the above eqnation dr is the dimension of the representation r and Xr{U) the 
character of a matrix U in this representation; the coefficients Dr{l3t) are given by 


DriP) 


FriP) 

drFoiP) ’ 


(13) 


with 

/ OO 

^ detIr^+,+niN/3) . (14) 

n=—oo 

The Tj’s are a set of integers labelling the representation r and they are constrained 
by ri > • • • > r^v = 0. The indices 1 < j < iV label the entries of the N x N 

matrix the determinant of which is taken and In{P) denotes the modified Bessel 

fnnction of order n. 

As a result of neglecting the space-like plaquettes, each space link variable ap¬ 
pears only in two time-like plaquettes, and the integration over the variables Ui{x, t) 
can be easily performed by using the rules for the integration of characters (see Q] 
for the details). One gets the following effective action for the Polyakov loops: 

e^P°'(^^=o) =nE[^-(A)]'V‘x.(P(f+i))x.(Pt(f)) . (15) 

x,i P 


From the effective action we obtain some general informations on the scaling be¬ 
haviour that we may expect for the critical coupling. They are obtained by requiring 
that the physics is unaffected by changes of Nt provided fdt is suitably rescaled. In 
particular it is easy to derive that the rescaling of f3t is linear. 

In fact if we just assume that fdt becomes large for large W, we can replace in 
(|I^) the coefficients Dr{(dt) with their large fdt behaviour which is well known and 
is given by 

lim Dr{(dt) = e , (16) 

Pt^oo 


where Cr denotes the quadratic Casimir in the representation: 


N N 

C'r = E rf - ^ - 

i=l i=l 

"^For the derivation of this and of the following formulas, see for instance Ref. 


V(iV2 - 1) 
12 


( 17 ) 
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So we have at the leading order 


\DA.I},)f‘ exp 


( CM \ 

1 2JVA(iV,)h 


(18) 


Clearly the above discussion is independent from the asymmetry of the lattice, and 
in particular the case of the symmetric lattice can be recovered by simply replacing 
/3t and Nf with (3 and rit respectively. The r.h.s. of ([T8|) , and hence the effective 
action (^), are independent from iV* in the large Nt limit provided Pt^Nt) depends 
linearly from Nt- If we dehne, consistently with eq. 


J 


lim 

Nt^oo 


mt) 

Nt 


(19) 


we obtain the large Nt limitQ of the effective action (^): 


Qhk 

e 


nEx r{P{x + i))xr(T’''’(x)) exp 

X,l ^ 


a 


2NJ 


( 20 ) 


The above effective action is the same one that we would have obtained, for any 
value of Nt and not just in the large Nt limit, if we had started, always within the 
f3s = 0 approximation, from the heat kernel action instead of the Wilson action and 
with a coupling /3hk(-X^t) given by 


MNt) = JNt . 


( 21 ) 


The critical value of J at the deconhnement transition will be determined in 
Sec. H by using an improved mean held method and a Poisson resummation of eq. 
(^). In the next section instead we shall calculate the next to leading term in 
the rescaling of (3t, namely the coefficient ai of eq. (||). Unlike J, this depends on 
the particular form of the action on the lattice (it is zero if we start from the heat 
kernel action as shown in (pT])) and it is completely determined by the form of the 
coefficients Dj.. 


3 The scaling law: exact calculation of the sub¬ 
leading term 

We have shown in the last section that the effective action ([T5| ) has a smooth large 
Nt limit, given in eq. (^0|) , provided a linear rescaling of f3t is assumed. In this 
section we want to go further, and assuming a linear riscaling of the form 

_ (3t{Nt) = JNt + ai (22) 

®The large Nt limit can correspond either to the continuum limit, for instance in the symmetric 
lattice where p = 1 and Nt = rit, or to the limit p ^ oo, which is a continuum limit only in the 
time dimension. These two limits, although conceptually completely different, produce here the 
same effect as a consequence of having neglected the space-like plaquettes. 
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we shall determine the next to leading coefficient ai by requiring that the effective 
action is independent from also in the next to leading order of its expansion. 
More precisely we shall show that it exists, for each N, a value of oi, independent 
from the representation r, for which is independent from Nt in the large 

Nt limit up to order 

In order to do that we need to find the large fdt asymptotic expansion of Dr{(3t) 
up to the order -4, namely one order higher than the leading one, already given in 


eq. (0) and used in the derivation of the heat kernel action (^). 
To begin with let us write Fr{j3), defined in (|^, as 


FriP) 







(23) 


The infinite sum over n implies the costraint = 0; as a consequence a shift of 
all the integers by an arbitrary, not necessarily integer, quantity does not affect 
the r.h.s. of (^31). We can therefore replace the rfs in (^) with 


N 




k=l 


which satisfy the relation 


= 0 • 


3 


We then proceed to the following redefinitions: 

fi 




^3 = 


' 3 


n 


V = 


With these replacements eq. (^) can be rewritten as 

F,{I3) = '£ n e vW , 

. J-^x/N3 . \/N i 


where is the unitary Vandermonde determinant defined by 

J(«) = n2smh^. 

i<j ^ 


(24) 

(26) 

(26) 

(27) 

(28) 


We shall now derive the asymptotic expansion in ^ of (^), keeping Vj fixed, in 
spite of the (3 dependence hidden in the definition of Vj. Only at the end of the 
calculation we shall replace Vj with its expression. 

The asymptotic expansion of the integral at the r.h.s. of (|^) is obtained in 
two steps. First, by noticing that under the integral = —i^, we extract the 
Vandermonde determinant and write 


Fr{fi) = J{ 


— 1 


d 


E 


T^/WP 


\/W[3 du V 




^iJ2fiFj+u)^j+NpJ2.cos^ 


(29) 
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The asymptotic behaviour of the above integral can now be easily calculated by 
expanding the cosine and using the saddle point method. Also, the integrals over the 
can be taken from —cxo to +cxo up to terms which are exponentially depressed. 
The result is: 


FrW) = 




X 


'y/Nf3 dv' 


where C collects all terms that do not depend on the representation, that is that 
are independent from Vj. These are in fact irrelevant for our calculation as only 
the ratio appears in (plSf ). Similarly the terms of order have not been 
evaluated because they are either representation independent or of order higher 
than dj once the (3 dependence of Vj is taken into account. 

The next step is the summation over v. It is easy to see that up to non pertur¬ 
bative terms we can perform the replacement 

^ ^ r du (30) 

V J-oo 


and that the resulting gaussian integrals can be calculated. With the condition 
Yjj = 0, the result is: 


FriP) = J{ 


-i d ^ iy ^2 

3 2 3 








(31) 

We have now to perform the derivatives included in the unitary Vandermonde de¬ 
terminant i7. By expanding the sines (see eq. (p8|)) we can first write it, up to the 
relevant order in as 


Ji 


-i a. 




-i d 




), 


(32) 


where A is the ordinary Vandermonde determinant]^. By operating with the r.h.s. 
of (^ we get 


FM = A( 


-i 9 , _iy ,^.2 

2 3 




)e' 


1+ i:d+377Fi;i:d+g 


24/3 A2 


3 


24:Np 


(33) 


®In deriving eq. (^) the relation J2j ^ which follows from J2j been taken 

into account. 
































Finally we have to operate with A( 


-i d 


) and keep in the result only the terms 


quadratic in vj. In fact terms independent from Vj are irrelevant for the reason 
given above, and terms of the type Yi,j are of order once the dependence 
of Uj from [3 is taken into account. We shall make use of the following relations: 


A(- 

A(- 


A( 

-) 

dp’ 

A) 

dp’ 


= A(0e-4E 


dp 


2 

J "3 




.Z/2 

3 3 


AW)('£,v. 




^ 4 


3 


A(0(i: Pj - MN --)Y. v] + C(iV))e-3E.-] 

0 j 


(34) 

(36) 

2 

; 

(36) 


Eq. (^) simply follows from the remark that the polynomial in front of the expo¬ 
nential is globally of order jg completely antisymmetric under 

exchange of the p/s. Eq. (|35| ) is a simple consequence of (|3^): rescale the Pj^s 
in the exponent by a parameter 7 and take on both sides the derivative respect 
to 7. Eq. (p 6 D can be checked directly for small values of N and then proved by 
induction. Notice that in deriving eq. (|36|) the relation J2j = 0 has been used. 
By applying the above relations to eq. ( P3[ ) we hnally obtain 

FriP) = . (37) 


It is now sufficient to insert (^) into (pAf ) and to keep into account the dehnition of 
Pj given in and the expression (|I^) of the quadratic Casimir in order to obtain 


DPP) = e 


Cr 
■ 2N 


1 I N^-2 1 

0^ N'^ 


+0(/3-3) 


(38) 


This asymptotic behaviour can be inserted into the effective action with the 
rescaling of the coupling given in eq. (p^) to determine its large Nt behaviour: 


g5poi(/3s— 0 ) 


nEe 



Ar ^-2 

4JV^ 



Xr(-P(S+i)),\r(f’’(S)) . 


(39) 


The leading order in gives, as we already knew, the heat-kernel action. As 
discussed earlier on, the rescaling of the coupling constant is determined by the 
requirement that the r.h.s. of (|3^) is Arindependent up to the terms O(^). This 

hxes the value of oi; in fact the term in the coefficients of the character expansion 
vanishes, for any representation, provided 

N^-2 
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Let us remark that the above result in the large N limit (ui = 1/4) was already 
obtained in ||^ by using an Eguchi-Kawai reduction scheme. 

We have determined, in a purely analytic fashion, one of the two coefficients 
appearing in the scaling law In the symmetric lattice (W = rit) case eq. 

(p^) describes the scaling towards the continuum limit (n* —^ cxo) of the critical 
temperature. This is the prediction that we are most interested in, since it can 
be tested againts the Monte Carlo data, all of which are performed on symmetric 
lattices. 

For asymmetric lattices (W = prit) eq. gives the relation between the 

couplings at different values of the asymmetry parameter^ p\ 

(3t{pnt) = Jpnt + ai = p - ai) + ai . (41) 

The last equation relates the value of Oi which we just obtained with the parameters 
Q,o,i introduced in (|T^), and it gives 

al = ai , «r = • (42) 

It is interesting to compare our prediction with the results obtained by Karsch 
0 in 3 + 1 dimensions and including space-like plaquettes. Using the same normal¬ 
ization, they are given by 

«+-( 0 . 1305 -^) = -01 + 0 . 1195 . ( 43 ) 

Our prediction for the subleading term a/ coincides with Karsch’s results, while 
a correction, presumably due to contribution of the space-like plaquettes appears 
in a^. 


4 The leading term in the scaling law and com¬ 
parison with the Monte Carlo data 


Unlike the subleading term ai calculated in the previous section, the critical value 
Jc of the rescaled coupling J depends on the dynamics of the system and cannot be 
obtained exactly, since it would require solving the model dehned by eq. An 

estimate of Jc can be obtained by using a mean held approximation. With the form 
of the effective action given in eq. ( ]20| ) this is not an easy task, due to the inhnite 
sum over the group representations appearing in eq. (^). However the sum over 
the representations can be done explicitely by means of a Poisson resummation, 
and the action can be recast in the following form (see 1,0): 




n + 


x,{i=l,2} 


47r/ 


exp 


W2-U 

24J 


E 


(iVJ)d 


-N)/2 


j{e{x))j{e{x + i)) 


(44) 


^Eq. (g| has the same form as eq. 
plaquettes, the quantum corrections — oi,ai in 


however, since we are neglecting the space-like 
dH) can differ from , a\ in (©■ 
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X 


{li} 


NT ^ / 
2 = 1 


x) 


■9 


'p(i)(f+i) +27r/i) 


where 9i{x) are the invariant angles of the Polyakov lines: P{x) = diag{e‘^^^^}; 
JT (6) is the Vandermonde determinant for a unitary matrix dehned in (P^D and P 
denotes a permutation of the indices. The set of integers {k} are winding numbers, 
related to the fact that the angles 6i{x) are defined modulo 27r. For special unitary 
groups the sum of the invariant angles must vanish modulo 27r; if we choose them 
to satisfy the relation + i) = 0, then the relation J2ik = ^ also 

follows. 

An inhnite sum over the integers {/*} also appears in eq. (|4^) , but the con¬ 
tribution of the winding modes is known to vanish in the deconhned phase in the 
large N limit and it is exponentially suppressed for finite N. Hence an excellent 
approximation is obtained by truncating the sum keeping only very low values of 
{li}. The mean field analysis is greatly simplihed in comparison with the original 
form (pOf) of the effective action, to the extent that improvements of the mean field 
approximation, obtained by considering larger and larger clusters of “spins” (see 
Fig. 1), become viable. In this way more precise estimates of the critical coupling 
can be obtained (see [Q for further details). The results of our mean field analysis 
for the SU(2) and SU(3) models are reported in Tab. I together with the N = oo 
result already obtained in |^. Notice that the complexity of the analysis increases 
with N] so we could reach a “step 3” cluster in the SU(2) case, while we could only 
go as far as a “step 2” cluster in the SU(3) case. Notice however that as the num¬ 
ber of internal degrees of freedom increases the mean field approximation becomes 
more and more precise (in the large N limit it is exact), and smaller clusters are 
sufficient to extract rather good approximations of the exact result. In the SU(2) 
case, having at our disposal three steps we tried a (naive) extrapolation of our es¬ 
timate toward the infinite cluster limit, assuming a simple geometric progression. 
This result is reported in the fourth line of Tab. I, while in the last line we have 
reported, for comparison, the Monte Carlo estimates of Jc (that we shall discuss 
below) for iV = 2 and N = 3. 


4.1 Comparison with the Monte Carlo results 

It is interesting to compare our predictions for Jc and ai with the existing Monte 
Carlo data in (2-1-1) dimensions. Very precise estimates of the deconhnement tem¬ 
perature exist for the N = 2 and N = 3 models in the range 2 < nt < Q. These 
can be found in |P, [1^, HI and are reported in Tab. II. All these simulations were 


made with the standard Wilson action and at p = 1. 

The data show that the expected linear dependence from nt is very well fulfilled, 
and, what is more important, they are precise enough to allow to measure, besides 
the leading linear dependence, also the subleading correction. Both in the N = 2 
and in the N = 3 case the Monte Carlo data fit the linear scaling law (eq.(^) for 
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Table I: Results for Jc obtained by (improved) mean field techniques For the N = oo 
case, see The boldface entries represent our best estimates of the value of Jc- 
In the last line we have reported the Monte Carlo estimates of Jc for N = 2 and 
N = 3. 



N = 2 

N = 3 N = oo 

step 1 

0.292 

0.332 0.351 

step 2 

0.332 

0.359 

step 3 

0.355 


extrapol. 

0.367 


Monte Carlo 

0.380(3) 

0.366(2) 


Table II: The critical coupling Pc as a function of the lattice size in the t direction, 
nt, in the (2+1) dimensional SU(2) and SU(3) LGT, taken from JI1|/ and 


nt 

N = 2 

iV = 3 

2 

~ 0.866 

0.906(2) 

3 

~ 1.251 


4 

1.630(8) 

1.638(6) 

6 

2.388(10) 

2.371(17) 


symmetric lattices) 


^ JcUt + a. 


(46) 


with a good confidence level. If we extract from the data the best fit valnes of 
Oi and Jc we obtain the resnlts reported in Table III. The comparison with onr 
theoretical prediction shows a qnite remarkable agreement even in the subleading 
term. 

Another way to compare the results of the Monte Carlo simulations with our 
analytic estimates is to factor out the leading linear dependence by plotting 
for various values of nt- The results are shown, for both Monte Carlo and analytic 
data, in Table IV and plotted in Fig. 2. The analytic data in Table IV are obtained 
by inserting in (^Sj) our best estimate of Jc and the exact result for Oi. The last 


row in Table IV is obtained from the best fit of the Monte Carlo data, using 
At the end of Sec. 1 we formulated the conjecture that in 2 + 1 dimensional LGT 
the dynamics of the Polyakov loop is dominated by the time-like plaquettes, and 
that the space-like plaquettes can therefore be neglected. This conjecture was based 
solely on the fact that the time-like plaquettes dominate for very low values of rit 
and that the scaling law for the critical coupling is the same for the theory with and 
without space-like plaquettes. The results of our analysis show that this conjecture 
al least predicts with very good approximation the location of the deconfinement 
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Table III: Best fit values of Jc and ai compared to our analytic predictions. 


N=2 N=3 

Fit of M.C. data Analytic Fit of M.C. data Analytic 

Jc 0.380(3) 0.367 0.366(2) 0.359 

ai 0.106(1) 0.125 0.174(6) 0.194 

Table IV: MC data for f3c{nt)/nt compared to our analytic predictions obtained 
inserting our best estimates for Jc and our exact result ai in eq. (i- Entries in the 
last row are obtained by fitting the MC data (see Tab. Ill) 



N = 

2 

N = 

3 

N = oo 

nt 

M.C. data 

Analytic 

M.C. data 

Analytic 

Analytic 

2 

~ 0.433 

0.430 

0.4531(8) 

0.456 

0.476 

3 

~ 0.417 

0.409 


0.424 

0.434 

4 

0.4075(19) 

0.398 

0.4094(14) 

0.408 

0.413 

6 

0.3979(16) 

0.388 

0.3952(27) 

0.391 

0.393 

oo 

0.380(3) 

0.367 

0.366(2) 

0.359 

0.351 


transition. The problem of whether the order of the phase transition and its critical 
exponents can also be correctly predicted in the same approximation is an open one 
and it wonld deserve fnrther investigation. 
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Step 3 


Fig. 1: Progressively more precise Bethe-like approximations. In Step 1 (standard 
mean field approximation) we integrate only on the invariant angles 6i of the central 
site, while all the other are kept fixed to the mean value 9 and the consistency 
condition is {P{9)) = P{9), where P is the Polyakov loop. In the successive steps 
also the invariant angles 9[ of the nearest neibourghs ("Step 2 ) and 0" of the next to 
nearest neibourghs ("Step 3) are integrated over, and the consistency condition reads 
{P{9)) = {P{6')). A link factor of the action (^) is associated to each link drawn 
with a full line in the figure. 



Fig. 2 : The avaliabe M. C data for f5c{nt)/n{t) are compared to scaling laws obtained 
in the present paper. 
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